Detection of supernovae neutrinos with neutrino-iron scattering 
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The f e - 56 Fe cross section is evaluated in the projected quasiparticle random phase approximation 
(PQRPA). This model solves the puzzle observed in RPA for nuclei with mass around 12 C, because 
it is the only RPA model that treats the Pauli principle correctly. The cross sections as a function of 
the incident neutrino energy are compared with recent theoretical calculations of similar models. The 
average cross section weighted with the flux spectrum yields a good agreement with the experimental 
data. The expected number of events in the detection of supernova neutrinos is calculated for the 
LVD detector leading to an upper limit for the electron neutrino energy of particular importance in 
this experiment. 

PACS numbers: 23.40.-s, 25.30.Pt, 26.50,+x 



A careful knowledge of the semileptonic weak inter- 
actions in nuclei allows the possibility of testing impli- 
cations of physics beyond the standard model, such as 
exotic properties of neutrino oscillations and massive- 
ness. The dynamics of supernova collapse and explo- 
sions as well as the synthesis of heavy nuclei are strongly 
dominated by neutrinos. For example, neutrinos carry 
away about 99% of gravitational binding energy in the 
core-collapse of a massive star, and only a small fraction 
(~ 1%) is transferred to the stalled shock front, creat- 
ing ejected neutrino fluxes observed in supernova rem- 
nants 

It was shown in Ref. [2| that accurate nuclear struc- 
ture calculations are essential to constrain the neutrino 
oscillations parameters of the LSND experiment [1] . This 
was also noted in previous works, e.g. Hayes et al. in 
Ref. 0]. In that work, based on a shell model, the same 
exclusive cross section in 12 C is obtained as with another 
shell model used in Ref @. This shows the importance 
of including configuration mixing (as done in both ref- 
erences) for this nucleus [37|- Nevertheless the QRPA 
predictions of Ref. [5j do not yield good results for this 
nucleus because the configuration mixing is not properly 
accounted for and the projection procedure (as done in 
Ref. 0) is not included. In particular, the employment of 
PQRPA for the inclusive l2 C(v e , e~) 12 N cross section, in- 
stead of the continuum RPA (CRPA) used by the LSND 



collaboration in the analysis of za, 



oscillations of the 



1993-1995 data sample, leads to an increased oscillation 
probability. Then, the previously found consistency be- 
tween the (sin 2 28, Am 2 ) confidence level regions for the 
Vjx —* Ve and the — > v e oscillations is decreased [2| . 

The measured observables are flux-averaged cross sec- 
tions. The KARMEN Collaboration measured charged 
and neutral cross sections induced on 12 C Q. They also 
measured (the only experimental data for a medium- 
heavy nucleus) the neutrino reaction 56 Fe(^ e , e~) 56 Co 
from e _ -bremsstrahlung with the detector surrounding 
shield 9]. This cross section is important to test the 
ability of nuclear models in explaining reactions on nu- 
clei with masses around iron, which play an important 



role in supernova collapse llOf . Experiments on neutrino 
oscillations such as MINOS [llj use iron as material de- 
tector, and future experiments, such as SNS at ORNL 
[l2| plan to use the same material. 

In a recent work, Agafonova et al. [l3j studied the ef- 
fect of neutrino oscillations on the supernova neutrino 
signal with the LVD detector. This detector studies su- 
pernova neutrinos through their interactions with pro- 
tons and carbon nuclei in a liquid scintillator and with 
iron nuclei in the support structure. Several estimates on 
deviations of the detected signal arising from different 
constraints on the astrophysical parameters, oscillation 
parameters and the non-thermal nature of the neutrino 
fluxes were studied before [l3[. Nevertheless, in all their 
estimates the corresponding i/-nucleus cross sections were 
kept within strict limits. 



In this work, we calculate the 



_56 



Fe cross sections 



using QRPA and PQRPA models to account for allowed 
and forbidden transitions. The present calculations are 
the first within the PQRPA framework for this purpose. 
In Ref. [IH it was established that PQRPA is the proper 
theory to treat both short range pairing and long range 
random-phase (RPA) correlations. When QRPA was im- 
plemented for the triad { 12 B, 12 C, 12 N} there were diffi- 
culties in choosing the ground state of 12 N because the 
lowest state is not the most collective @. PQRPA solves 
this puzzle because it treats correctly the Pauli Principle, 
yielding better results for the distribution of the Gamow- 
Teller (GT) strength. This problem does not exist in 
heavier nuclei, where the neutron excess allows QRPA to 
account for pairing and RPA correlations [l5[ . In the case 
of medium-heavy nuclei, such as 56 Fe, the consequences 
of the projection technique procedure can be manifest. 

Many calculations of the 56 Fe(z^ e , e~) 56 Co cross sec- 
tions with microscopic and global models have been re- 
ported previously. The first were shell model calculations 
developed by Bugaev et al. fig . They obtained the v- 
nucleus cross sections as a function of the incident neu- 
trino energy. A second estimate was obtained by Kolbc 
et al. [T?! using a nuclear Hybrid model: shell model for 
the GT and Fermi (F) transitions, and continuum RPA 
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(CRPA) for forbidden transitions. This cross section was 
employed to estimate the number of events from i/- 56 Fe 
reactions in the LVD detector [l|. Lazauskas et al. [IH 
used QRPA with the Skyrme force to explore the possibil- 
ity of performing nuclear structure studies using neutri- 
nos from low energy beta-beams Several ^-nucleus 
cross sections for different nuclei were also obtained re- 
cently with the relativistic QRPA (RQRPA) [20]. The 
v e — 56 Fe cross sections were also described with the gross 
theory of beta decay (GTBD) [p. 

The cross section for v e + (Z, A) — > (Z + 1, A) + e~ is 
given by 

a(E e ,J f ) = \v&F{Z + l,E e ) J d{ C o S 6)%{\h\,J } ), 

(1) 

where F(Z + 1, E e ) is the usual scattering Fermi func- 
tion, k — p e — q v is the momentum transfer, p e and q v 
are the corresponding electron and neutrino momenta, 
and 6 = q„ • p e is the angle between the incident neutrino 
and emerging electron. The <r{E e , Jf) cross section is ob- 
tained within first-order perturbation theory according 
to Ref. [7], where velocity-dependent terms are included 
in the weak effective Hamiltonian. The transition ampli- 
tude 7^(|k|, Jf) depends on the neutrino leptonic traces 
and on the nuclear matrix elements (NME), as explained 
in Ref. 0- Here, the NME are evaluated in QRPA and 
in PQRPA. We employ the (5-interaction (in MeV fm 3 ) 

V = -4n(v s P s + v t P t )5(r), 

with different coupling constants v s and Vt for the 

particle-hole, particle-particle, and pairing channels. 

This interaction was used in Refs. [22L I23L I24I |25| leading 
to a good description of single and double /3-decays. 

For 56 Fe we work within a configuration space of 12 
single-particle levels, including the oscillator shells 2hoj, 
3hu and AJvjj. The single-particle energies of the active 
3Hu> shell correspond to the experimental energies of 56 Ni. 
For the other 2fkv and Ahw shells we have used the har- 
monic oscillator energies with Hui/MeV = 45 A 1 / 3 — 
25 A 2 / 3 . The parameters v pair (p) and v pair (n) were ob- 
tained with the procedure of Ref. [26], i.e., by fitting the 
experimental gap pairing energies of protons and neu- 
trons, A n>p (N,Z) (eq.(2.96) of Ref. [13), to A„, p defined 
by the usual BCS equations. The BCS or PBCS equa- 
tions were solved in the full space of three oscillator shells. 
For the particle-hole channel we have used v ph = 27 and 
vf h = 64 (in MeV fm 3 ). These values were fitted to 48 Ca 
from a systematic study of the GT resonances [25| and 
shown to yield a good description of double beta decay. 
For the particle-particle channel, it is convenient to define 
the parameters 

_ 2yPP 

s ~ vr r ( P )+vr r (n)' 



and 

2vf 

pair 1 \ . pair / \ ' 
Vs (P) + v% (n) 

associated to the coupling constant of the T = 1,5 = 
(singlet) and T = 0, S = 1 (triplet) channels respectively 
Ref. [2a]. We adopt s ~ 1, which restores the isospin 
symmetry in QRPA for N > Z. As the experimental 
errors in the averaged cross sections are very large, the 
agreement of the theoretical cross section is not suffi- 
cient to select the best nuclear structure calculation and 
other observables must be found. We use the behavior 
of the B(GT~) strength as function of the parameter t 
to conclude that better results could be obtained when 
the particle-particle channel is off, t = 0. With this value 
the theoretical value(i?(GT _ ) = 17.7) overestimates the 
experimental value (9.9 ± 2.4 [2§]) similarly to previ- 
ous and more sophisticated QRPA calculations for 56 Fe 
(B(GT-) = 18.68 [13) with the Skyrme force. 
The flux averaged inclusive cross section reads 

(<i e ) = J dE v a{E v )n{E v ), (2) 

where a e (E v ) — J2jj a e(E e ,Jf), is the inclusive cross 
section as a function of the neutrino energy and n(E v ) is 
the neutrino normalized flux. As a first test, we fold the 
<j e {E v ) with the Michel energy spectrum 

n{E v ) = -^{M li -2E v ), (3) 

where Af M is the muon mass. In Table U we compare our 
56 Fe(z/ e , e") 56 Co cross section (tr e ) in QRPA and PQRPA 
with other nuclear models for the energy window of fi- 
Decay-At-Rest (DAR) neutrinos that the KARMEN ex- 
periment observed. 

TABLE I: Comparison of (<j e ) in 10" 42 cm 2 for 
56 Fe(^e,e-) 56 Co obtained in our QRPA and PQRPA 
confronted with other nuclear models. For the Hybrid 
model [ItJ we denote with (a) partial occupation in the 
ground state and with (b) no occupation. 



Model 


{^) 


QRPA 


264.6 


PQRPA 


197.3 


Hybrid (a) [17] 


228.9 


Hybrid (6) [17] 


238.1 


TM [31J 


214 


RPA [32] 


277 


QRPA S [18] 


352 


RQRPA [20] 


140 


Exp [9] 


256 ± 108 ± 43 



From Table[T]we note that our results for (<r e ) — 264.6 x 
10~ 42 cm 2 (QRPA), and 197.3 x 10~ 42 cm 2 (PQRPA) are 
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FIG. 1: (Color online) Inclusive 56 Fe(> e , e") 56 Co cross sec- 
tions (in 10 -42 cm 2 ) evaluated in QRPA (dashed line) and 
PQRPA (solid line), in the DAR region, are compared with 
those obtained with other nuclear structure models: GTBD 
(dashed-dot-dot-dot line) [HI, Hybrid (dashed-dot line) [2S| |. 
QRPA S (dashed-dot dot line) [II], and RQRPA(dashed 
line) SI. 




40 60 

E v (MeV) 

FIG. 2: (Color online) Same as in Figure [1] The inclusive 
cross sections (in 10 -42 cm 2 ) are shown. In this case, the 
neutrino energy window is characteristic of the LVD experi- 
ment. 



flavor a, it is 



in agreement with the experimental value 256±83(siai)± 
42(syst) x 10~ 42 cm 2 Q. The main difference between 
QRPA and PQRPA, both solved consistently with the 
same interaction, shows that the projection procedure is 
important in a medium mass nucleus such as 56 Fe. 

In FigureQ]we plot the inclusive 56 Fe(^ e , e~) 56 Co cross 
sections (in 10 -42 cm 2 ) evaluated in QRPA (dashed line) 
and PQRPA (solid line), in the DAR region. A compar- 
ison is shown with other nuclear structure models. All 
models yield the same energy dependence that goes ap- 
proximately as £" 2 for low incident neutrino energies, ex- 
cept for the GTBD model, which shows a large deviation 
from the other cross sections. 

Figure [2] excludes the GTBD results and extends the 
energy scale to 100 MeV (supernova neutrino energies). 
a(E v ) reaches the DAR energy region at E v ~ 60 MeV. 
Beyond that the QRPA result is above the other models. 
Nevertheless, the PQRPA cross section lies below, show- 
ing the effect of the projection procedure. In this region, 
the main contribution arises from the non-allowed tran- 
sitions, as found in previous works [l8ll20f. 

The number of events detected for supernova explo- 
sions is calculated as, 



L a 



El 



N a =N t T a (E v ) ■ a(E v ) ■ e{E y )dE v 



(4) 



where the index a = v e ,u e ,v x and (y x = i/ Tl v tll v fJj ,i' T ) 
indicates the neutrino or antineutrino type, N t is the 
number of target nuclei, J- a (E u ) is the neutrino flux, 
a(E v ) is the neutrino- nucleus cross section, t{E v ) is the 
detection efficiency, and E v is the neutrino energy. Re- 
cent calculations by the LVD group [l3[ estimate that 
the [v e + V e ) interactions on 56 Fe are almost 17% of the 
total detected signal. 

The time-spectra can be approximated by the zero- 
pinched Fermi-Dirac distribution. For the neutrino of 



4ttD 2 T£F 3 (0) exp(E v /T a ) 



1 



(5) 



where D is the distance to the supernova, E u is the 
neutrino energy, L a is the time-integrated energy of fla- 
vor is ai T a is the neutrino effective temperature, and 
^3(0) = / °° d 3 x x 3 /(e x + 1) is the normalization factor. 
For a galactic supernova explosion at a typical distance 
D = 10 kpc, it was assumed that the total binding en- 
ergy for each flavor is L a = f„ a Eb, with Eb = 3 x 10 53 
erg, and a perfect energy equipartition between the neu- 
trino flavors, f a = fv e = fp e = fv x = 1/6. Hence, it 
is possible to assume that the fluxes ( v ^ , v T , v ^ , v T ) are 
identical. Because the pinched factor was assumed zero 
for all neutrino flavors, we can fix the effective neutrino 
temperature as T Vm — 1.5T &e and T„ e = 0.8T Pe , leaving 
Tp e as a variable parameter in the interval studied in 
Ref. pj]. 

When the neutrinos escape from the star, they cross re- 
gions of different densities where a flavor transition could 
happen. Usually one assumes two resonance layers which 
we call Mikheyev-Smirnov-Wolfenstein (MSW) reso- 
nances throughout the text (see for example Ref. [33|). 
According to the mass scheme shown in [l3| , the observed 
electron neutrino fluxes (T v <, , ) originating from MSW 
resonances are linear combinations of the original neu- 
trinos fluxes in the star, J 7 ^ and , with coefficients 
governed by the crossing probability in the high density 
resonance layer, P fJ (A 2 tm , #13). For simplicity, we only 
show differences that appear in the number of events cal- 
culated from the convolution of cross sections obtained 
with different nuclear structure models with the original 
supernova fluxes, i.e., 

/>oo 

N e ee N e {T v J = N t / T Q e {E Vl T„ c ) ■ a e {E v ) ■ e{E v )dE v , 
Jo 

(6) 
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for "direct" electron neutrino event, and 



N e = N e (T Vx ) = N t 



/ ^°(E u ,T Ux )-a e (E v )-e(E u )dE u , 
Jo 

for the "indirect" number of events for electron neutrino 
associated to the total z/ e -flux coming from the contri- 
bution of J^.. Due to the MSW effect, electron neu- 
trino fluxes mix with non-electron neutrino fluxes (i.e., 
v x = v^^Vr), and therefore with the MSW resonance the 
Ve§ might get a "hot" contribution to their flux. An- 
other important issue, not considered for simplicity in 
the present work, is the spectral swapping of the neu- 
trino flux (Ref. [3J]). Duan et al. have shown that 
certain numerical results in the simulation of neutrino 
and antineutrino flavor evolution in the region above the 
post-supernova explosion proto-neutron star cannot be 
easily explained with the conventional MSW mechanism 
Ref. Q. 

For the neutrino reactions 56 Fe(z/ e , e _ ) 56 Co, we calcu- 
late N e and N e as a function of the neutrino tempera- 
tures T„ e and T Vm , folding a e (E u ) from different nuclear 
structure models with the neutrino fluxes (E v , T Ue ) 
and J-JJ (E u ,T Val ), respectively. The limits for the tem- 
peratures, T Ve and T Ux , were obtained from the interval 
T Pc G [4,7] MeV and the relations T Vtc = 1.5T Pe and 
T Vc = 0.8T Pc , employed by the LVD group [Hj|. The 
e(E u ) efficiency is taken from Figure fb of Ref. [HI]. 
The results obtained are shown in Figure [31 The left 
panel shows the number of events for electron neutri- 
nos, N e , with different ajE v ), our QRPA and PQRPA, 
QRPA S [jj], RQRPA and the Hybrid model [H 
employed by the LVD detector. On the right panel we 
show the number of events N e . Although one knows that 
v x neutrinos at supernova energies can only induce neu- 
tral current events, we evaluate this quantity because it 
will be modified by MSW oscillations according to the 
scheme presented in equations (2) and (4) of Ref. [f3j], or 
in equations (fO) and (f 2) of Ref. [35[. Despite it is cer- 
tain that the N e could be obtained from the expression 
for N e extending the interval of T Ve to cover the interval 
T Vx , this region (T„ c G [6, f0.5] MeV) of temperature for 
v e is not in agreement with the physical range depending 
on the neutrino transport that it is predicted by different 
supernova modelers, as such in Ref. [36] . 

We note that N e and N e increase with the temper- 
atures T„ c and T Vm . The increase for each N e follows 
the increase of the different er e . The contribution of the 
neutrino flux, J-® , in N e is strongly concentrated in the 
region below 60 MeV. This is because: (i) the mean neu- 
trino energy (E Ve ) of the flux varies from 10.1 to 17.6 
MeV approximately 38] ; and (ii) the contribution of the 
product of a e with the flux tail is not important. Notice 
that the ordering of the N e in Figure [3] is the same as 
the ordering of a e shown in Figure [TJ For example, the 
crossing between iV e 's of our QRPA and Hybrid model 
at T„ e ~ 4.8 MeV originates from the crossing of the 
corresponding a e at E v ~ 32 MeV as Figure [1] shows. 



- - - - Hybrid 

QRPA 

QRPA, 




3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 10.5 

T Ve (MeV) T Vx (MeV) 



FIG. 3: (Color online) Number of events obtained from 
the convolution of the neutrino fluxes with the cross sec- 
tion obtained with different nuclear structure models: Hybrid 
(dashed-dot line) [H], QRPA (dashed line), QRPA S (dashed- 
dot dot line) [lj|, PQRPA (solid line), and RQRPA(dashed 
line) 



The above behavior also applies to N e , but they are 
shifted according to the shift that has with respect 
to . This means that the main contribution to N e 
comes from the convolution of J-~ with a e in the energy 
interval [18.9,33.1] MeV where the larger energy flux of 
v x is concentrated. The right panel of Figure [3] shows 
additional crossings at T Vx ~ 10.5 MeV which is a re- 
sult of the corresponding crossings of <r e (QRPA-QRPA s ) 
at E v - 56 MeV and er e (PQRPA- RQRPA) at E v - 60 
MeV, shown in Figure [TJ We conclude that the relevant 
energy interval for the v e — 56 Fe reaction is E v < 60 MeV 
for the astrophysical parameters adopted in LVD. 

In summary, we have employed the projected QRPA to 
calculate the 56 Fe(^ e , e~) 56 Co cross section. The calcu- 
lated cross section is compared with a QRPA calculation 
with the same interaction showing that the projection 
procedure is important for medium mass nuclei. The 
cross section is also compared with other RPA and Hy- 
brid models. The PQRPA yields smaller cross section 
than almost all RPA models with exception of relativis- 
ts QRPA [H for E v < 60 MeV. Above this energy and 
up to E v = 100 MeV, the PQRPA leads to the small- 
est cross section. Therefore, we feel that a more detailed 
study of allowed and forbidden transitions in the region 
below E v = 100 MeV is imperative, both experimentally 
and theoretically. In particular, the region with E v < 60 
MeV is the most important for the LVD detector [I3 |. 
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In a future work we plan to include the MSW effect in 
the same way as was as done by Agafonova et al. and 
an explicit account of the uncertainties in the supernova 
neutrino flux will be considered. 
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